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We present first results for lattice simulations, on a single volume, of the low-lying spectrum of an 
SU(3) Yang-Mills gauge theory with Nf = 10 light fermions in the fundamental representation. Fits to the 
fermion mass dependence of various observables are found to be globally consistent with the hypothesis 
that this theory is within or just outside the strongly-coupled edge of the conformal window, with mass 
anomalous dimension 7* consistent with 1 over the range of scales simulated. We stress that we cannot 
rule out the possibility of spontaneous chiral-symmetry breaking at scales well below our infrared cutoff. 
We discuss important systematic effects, including finite-volume corrections, and consider directions for 
future improvement. 



PACS numbers: ll.lO.Hi, ll.15.Ha, 11.25.Hf, 12.60.Nz 

Introduction In three recent papers HHH, we studied 
the properties of an SU(3) gauge theory with Nj mass- 
less Dirac fermions in the fundamental representation as 
Nf increases from 2 to 6. We noted that the Nf = 2 simu- 
lations are in good agreement with measured QCD values, 
and that the A'^^ = 6 results indicate substantial enhance- 
ment of the chiral condensate [IJ. We also observed that the 
spectrum becomes more parity doubled and the S parame- 
ter per electroweak doublet decreases as Nf increases from 
2 to 6 |;2J. Most recently, we examined vr — vr scattering 
in these theories ||3l, noting that the low-energy scattering 
length decreases when Nf increases from 2 to 6. 

The Nf = 6 theory is thought to exhibit these trends be- 
cause Nf has increased toward a critical value N'j at which 
the theory transitions from confinement and chiral sym- 
metry breaking to infrared conformal behavior. The pre- 
cise value of N'j remains uncertain, but all recent evidence 
points to confinement and chiral-symmetry breaking for an 
SU{2>) theory with Nf = 8 Ii4n7j, while most recent anal- 
yses indicate that an SU{?>) theory with Nf = 12 is con- 
formal in the infrared S |5l IS HTSll (although some studies 
have concluded that it is spontaneously broken lfT4l[T5l ). 

Here, we summarize results of numerical lattice simu- 
lations for an SU{2>) gauge theory with 10 light Dirac 
fermions in the fundamental representation. This the- 
ory was recently studied by Yamada et al lfT6UT8l . who 
computed the running gauge coupling in the Schrodinger- 
functional scheme, finding evidence that the theory is con- 
formal in the infrared, dominated by a relatively strong 
fixed point and with a large mass anomalous dimension. 

We examine the particle spectrum of the Nf = 10 



theory, including the masses and decay constants of the 
lowest-lying pseudoscalar, vector and axial-vector states, 
and the masses of the lowest-lying nucleon state and its 
parity partner. We also consider the chiral condensate. 
We compute these quantities on the lattice using a finite 
fermion mass m, and analyze the behavior as m is extrap- 
olated toward zero. 

Motivated by the running coupling results, we focus first 
on the infrared conformal hypothesis. We argue that the 
masses, decay constants, and chiral condensate scale with 
m over a range of m values in a manner consistent with 
mass-deformed conformal perturbation theory [[9l [T9H22l . 
indicating that the theory is either inside or just below the 
edge of the conformal window. We cannot rule out the pos- 
sibility that the theory will exhibit an intrinsic confinement 
scale and spontaneous chiral-symmetry breaking scale at 
lower values of m. As m is decreased, finite- volume ef- 
fects become more important, and these effects must be 
considered carefully, particularly since our current results 
are restricted to a single volume. 

Simulation details Simulations are performed using 
domain-wall fermions and the Iwasaki improved gauge ac- 
tion ||23l . The domain wall formulation suppresses the chi- 
ral symmetry breaking associated with fermion discretiza- 
tion, and preserves flavor symmetry at finite lattice spac- 
ing. Dimensionful quantities are given in lattice units, with 
implicit dependence on the lattice spacing a. The lattice 
volume is set to 32"^ x 64, with the length of the fifth di- 
mension Ls = 16 and the domain wall height uiq = 1.8. 
We choose j3 = Q/ = 1.95, which lies on the weak- 
coupling side of a bulk phase transition, and leads to a 
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vector- meson mass 0.24 < My < 0.34 for the range of 
fermion masses m used in our analysis. This is similar to 
our previous studies at Nf = 2 and Nf = 6 [1,2]. Each 
ensemble generated contains approximately 1200 gauge 
configurations. 

Input fermion masses ruf = 0.01 to 0.03 are included 
in the simulations. At finite lattice spacing, even with 
= 0, the chiral symmetry is not exact, with the vio- 
lation captured in a residual mass m^cs- The total fermion 
mass m is then m = m,/ + m^cs- For our simulations, 
rrires ^ 0.0017, so that mics <C m/ for all values of nif. 

To study the effects of thermalization and fixed topolog- 
ical charge, we have generated gauge field configurations 
from both disordered and ordered starts for most ensem- 
bles. We note that an m/ = 0.005 ensemble was also 
generated, but did not show signs of adequate thermaliza- 
tion over the number of configurations generated, and we 
will not discuss those results further here. 

Combination of ordered/disordered data The topo- 
logical charge Q is observed to evolve very slowly on all 
ensembles, so that disordered starts generally remain in 
a sector with large net topological charge, while ordered 
starts are essentially stuck in a sector with zero net topo- 
logical charge. On the subset of our ensembles where ex- 
trapolation of Q for disordered starts is possible, we find 
that these effects are able to explain most of the observed 
discrepancy between observable values on the two ensem- 
bles. In a future study currently in progress, we will include 
these topological-charge corrections explicitly by measur- 
ing the topological susceptibility on all ensembles. 

Since neither of our evolutions at a given mass point has 
sufficient topological tunneling, we combine the results us- 
ing the difference to estimate a systematic error. For an 
observable O, determined on a pair of ensembles with the 
same physical parameters but with different initial states 
(ordered and disordered), we assume that the samples of 
O are sufficiently large within a given ensemble that the 
central limit theorem applies. Each distribution is then 
Gaussian, with mean and standard error (/ii, cti), (/i2, '^2)- 
There remains an unknown bias of the mean computed 
within a given topological sector with respect to the true 
mean jl, resulting from a properly-weighted distribution 
over all sectors. 

We expect that the true mean lies somewhere between 
III and /i2. To obtain a conservative estimate of the true 
distribution, we take a uniform distribution of width 5 and 
center p, = (^Ui + /i2)/2 to describe our knowledge of 
the bias-corrected mean. Convolving this uniform distribu- 
tion with the Gaussian sample distributions, the combined 
mean is given by = /x, while the variance is equal to 
al = \{(tI + ai) + \5'^ . We determine the width 6 by 
assuming that the width as a fraction of the mean is con- 
stant for a given observable as a function of light fermion 
mass, taking the maximum fractional difference over the 
values observed as the best estimate. The resulting 6 val- 
ues range from 5-15% for most observables, so that the 5 
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FIG. 1. Edinburgh- style plot comparing ratios of physical ob- 
servables for simulations of the Nf — 2 (red), N f = 6 (blue), 
and Nf — 10 (green) theories. The points are joined together 
in decreasing order in the fermion mass, with the heaviest and 
lightest masses denoted by a triangle and box respectively. The 
combinations of observables plotted are Mj^/Fp vs. Mp/Fp, 
chosen to clearly illustrate finite-volume corrections which act 
with opposite sign on bound-state masses and Fp. 



contribution to the variance is always significant. 

Finite-volume effects Our simulations so far have used 
a single lattice spatial volume, L/a = 32. Finite-volume 
effects can in principle be incorporated through a con- 
trolled expansion in some function of ML, with M the 
relevant mass scale for a particular framework as in Ref. 
[9|. Here we will discuss finite-volume effects qualita- 
tively, and argue that such effects should be relatively small 
on a subset of our results. 

An estimate of the significance of finite-volume correc- 
tions can be made by comparing ratios of physical observ- 
ables. Edinburgh-style plots, which were used frequently 
as a diagnostic tool in the early literature on lattice QCD 
ll24l . can be particularly useful for theories other than QCD 
llT2l . In Fig. [1] we compare the ratios of observables 
^n/Fp and Mp/Fp. In order to better illustrate the ex- 
pected behavior, we include our own spectrum measure- 
ments for Nf = 2 and 6 H] 13 for comparison. 

The combinations Mj^/Fp and Mp/Fp are chosen to 
clearly show finite-volume corrections. Decreasing the 
lattice volume tends to increase Af^r and Mp while si- 
multaneously decreasing Fp relative to their values in the 
infinite- volume limit, driving the points up and to the right 
in the plot as m — )^ 0. On the other hand, the infinite- 
volume scaling behavior for a chirally broken theory has 
the points moving to the left as m — t- 0, as the pion mass 
scales to zero. For the Nf = 6 case, the latter trend 
changes into the former around mf = 0.010. 

For the Nf = 10 data combined using the method de- 
scribed above, no movement is seen within errors for nif > 
0.020, consistent with the possibility that Nf = 10 can be 
described by a mass-deformed conformal expansion. In 
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this expansion, all scales vanish as a common power of 
m at leading order, so that ratios of observables remain 
fixed as m is varied. The ruf = 0.015 and mj = 0.010 
points show small and large displacements, respectively, 
which may indicate finite -volume effects becoming large 
atnif = 0.010. The fits described in the next section, for 
nif > 0.010, > 0.015 and > 0.020, are consistent with 
these conclusions. 

The lack of scaling with m of the points on this plot 
is not itself sufficient evidence that the mass-deformed 
conformal framework is an adequate description of the 
Nf = 10 results. For a theory in a heavy-quark regime 
with all observables showing linear dependence on m, the 
ratios shown are also expected to approach a constant value 
for large m. However, the consistency of this plot with ex- 
pectations from the conformal hypothesis motivates further 
study within that framework. It is also clear from Fig. [T] 
that chiral perturbation theory cannot provide an accurate 
description of our Nf = 10 results: there are no signs of 
the expected decrease of Mp with respect to Fp, and the 
numerical size of Mp /Fp is too large for the chiral expan- 
sion to be convergent. 

Infrared conformal hypothesis If the Nf = 10 theory 
is conformal in the infrared limit, then the chiral symmetry 
is broken only by the explicit fermion mass, and the ex- 
pected mass dependence of the spectrum is determined by 
the emergent conformal symmetry. This possibility is sup- 
ported by evidence from running-coupling studies llT6l[T7l 
that the gauge coupling g^(/u) evolves slowly at long dis- 
tances. Therefore, we fit the spectrum using the infrared 
conformal hypothesis, assuming that g^{^) remains at its 
fixed-point value g*^, and that the data can therefore be de- 
scribed by mass-deformed conformal field theory ||9l [T9l - 
l22l with mass anomalous dimension 7*. 

When an explicit fermion mass m = m(A) ^ A is 
introduced, the running mass for scales below A is given by 
m{ij) = m{A/fi)'' . For some energy scale M <C A ~ 
1/a, the running mass thus satisfies the equality m{M) = 
M, so that at scales small compared to M the fermions 
decouple from the theory, leaving an effective pure-gauge 
theory which confines as the gauge coupling flows away 
from the fixed-point value g*"^. So long as g*"^ is reasonably 
strong, the induced confinement scale will be of order M. 
The mass of all fermion bound states is then given by ||9l 

Mx = Cx m[i/(i+T*)l + Dx m, (1) 

where we have included a small correction term. With the 
masses expressed in units of the cutoff A, Cx and Dx are 
dimensionless coefficients. Since the explicit breaking of 
chiral symmetry is of order M, there is no approximate chi- 
ral symmetry to be broken spontaneously. Thus this scaling 
law applies as well to the pseudoscalar mass. 

The pseudoscalar, vector, and axial-vector decay con- 
stants as we define them are expected to scale in the same 
way as bound-state masses \,22i . The chiral condensate has 
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FIG. 2. scans as a function of 7* for our Nf = 2 results (top) 
andA'^y — 10 results in the range mj > 0.015 (bottom). Contours 
shown are, from bottom to top at 7* = and Nf = 10: {ipTp) 
(red), AIx* (orange, dashed), Fp (purple), Af^ (blue, dashed), 
Fv (pink), Fa (pink, dashed), Mjy (orange). My (blue), Mp 
(green, dashed), and total (black). At Nf = 2, the pseu- 
doscalar mass shows the expected scaling behavior Mp ~ m, 
which appears as a minimum at 7* ^ 1 in this analysis. 

a more complicated dependence on the fermion mass ||9l: 

(^^/;) = Acm + Bcm^^^-^*^'^^+''^^^ 

+ Ccm[=^/(^+T*)l + (2) 

The above expressions vanish as m — t- with the scal- 
ing determined by a single parameter 7*, a behavior qual- 
itatively different from that of a theory with spontaneous 
chiral symmetry breaking. 

Results for global fits to the combined simulation data 
for the ranges > 0.010, mf > 0.015 and m/ > 0.020 
are shown in Tablejl] Since we have a relatively small num- 
ber of mass points to work with, we here consider only fits 
with the D-terms set to zero. As anticipated, the fit quality 
is reasonably good for the restrictions nif > 0.015 and 
ruf > 0.020. Including the = 0.010 data changes the 
fit parameters significantly, matching our expectation that 
finite-volume corrections become large for these points. 
The other two fits indicate a large anomalous dimension 
7* > 0.8 and are consistent with 7* = 1, the value antic- 
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Ubs. 


my > 0.010 


TTlf^ 0.015 


rrif > 0.020 


-A- 

7 


1.69(16) 


1.10(17) 


1.35(47) 


[68% CI] 


[1.54,1.86] 


[0.95,1.27] 


[1.06,1.73] 


[95% CI] 


[1.40,2.06] 


[0.82,1.46] 


[0.83,2.27] 


Cp 


0.98(9) 


1.44(21) 


1.21(37) 


Cv 


1.17(10) 


1.70(25) 


1.42(44) 


Ca 


1.43(13) 


2.14(32) 


1.79(56) 


Cjv 


1.75(16) 


2.53(37) 


2.10(65) 


Cn* 


2.23(25) 


3.35(55) 


2.87(92) 


Cpp 


0.121(12) 


0.190(28) 


r\ -I A / C -I \ 

0.164(51) 


Cpv 


0.165(15) 


0.238(35) 


0.195(60) 


CpA 


0.136(13) 


0.192(28) 


0.154(48) 


xVd.o.f. 


69/31 


14/23 


3.1/15 



TABLE I. Global fit results for the conformal hypothesis of 
Eqs. (|T]| and l[2]l, based on combined ordered/disordered data as 
described in the text. The labels P, V, A, N, N* correspond to the 
pseudoscalar, vector, axial-vector, nucleon and nucleon-prime, 
respectively. Decay constants for channel X are denoted by FX. 
Errors shown on all quantities are purely statistical, and ignore 
correlations between observables. For 7*, two-sided 68% and 
95% confidence intervals are also shown. The rrif > 0.010 fit 
(left column) has significantly worse /d.o.f., possibly due to 
the presence of finite-volume effects. 



ipated for a theory with Nf at the edge of the conformal 
window Il25ti27l . 

To better understand our fit results, we show in Fig. [2] 
scans over as a function of 7*, broken up for each indi- 
vidual observable included in the ruf > 0.015 fit. Several 
of the observables show an individual minimum in com- 
patible with the global best-fit value 7* ~ 1.10. The chiral 
condensate (shown in red) has no clear minimum, but it 
contributes very little to overall x^, so we omit it from the 
global fits in Table [l] The global fit with the condensate in- 
cluded is not significantly different, but exhibits very strong 
correlations between the parameters 7*, Ac and Be', this 
behavior is expected for 7"^ near 1 , for which the Ac and 
Be terms in Eq. (|2| are nearly degenerate. 

As is evident from Fig. |2] the distribution of x^/d.o.f. 
as a function of 7* is not symmetric about the minimum. 
We estimate a two-sided 68% (95%) confidence interval 
on 7* directly, varying by Ax^ = 1 (A^^ = 4) about the 
minimum of the x^ contour shown in Fig. [2] Results for 
each mass range are shown in Table |lj In all cases we find 
7* > 0.8 at two sigma. 

A similar plot using our Nf = 2 results is shown for 
comparison. As expected the Nf = 2 theory shows gen- 
erally very poor power-law fits for any 7* < 2, with the 
exception of the pseudoscalar mass (green, dashed), which 
scales as Mp ~ m in accordance with chiral perturbation 
theory. 
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FIG. 3. Simulation results for the pseudoscalar mass (circles), 
vector mass (squares), and axial-vector mass (triangles). Error 
bars on the points are estimated using the combination method 
described in the text. Two fit types are compared: linear My ~ 
ay + bym (red), and power law Aly ~ m^^^ (blue). The power- 
law fits correspond to a mass-deformed conformal theory with 
7* = 1. Only observable values with ruf > 0.015 (filled sym- 
bols) are used in the fits. 



Chirally Broken Hypothesis Despite the quality of fits 
obtained under the infrared-conformal hypothesis, it re- 
mains possible that the Nf = 10 theory is chirally broken. 
A rigorous test of this possibility would involve chiral per- 
turbation theory to extrapolate to m = 0. But as discussed 
in the context of Fig. [T] we do not expect this expansion 
to be convergent for my^ > 0.015. We have nevertheless 
attempted to fit our Nf = 10 results using NLO chiral per- 
turbation theory, as done previously for Nf = 2 ([Q and 
Nf = 6 li28l , finding (at Nf = 10) generally large values 
of x^/d.o.f. and best-fit values pointing to a poorly conver- 
gent expansion. We omit the details of these fits here, but 
will present them in a future work. 

An alternative, crude approach is to use the extrapola- 
tion formula Mp ~ bpm}/^ for the pseudoscalar mass, 
and the linear expression Mx ^ ax + bxm for the other 
masses and decay constants. In Fig. [3] we compare fits of 
this type for the vector and axial-vector masses to mass- 
deformed conformal fits with fixed 7* = 1, a value within 
the errors of our global conformal fit. Within the range of 
fermion masses considered, we cannot clearly distinguish 
this simple linear dependence from the power-law behavior 
of the mass-deformed conformal fits based on our current 
results. We also show fit results for the pseudoscalar mass; 
under either the conformal or chirally broken hypothesis, 
this state scales as A'lp ~ m^^^, so only a single fit is 
shown. 

Discussion We have presented here the first non- 
perturbative calculation of the spectrum for an SU(3) 
gauge theory with Nf = 10 light fermions in the fun- 
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damental representation. Previous studies of the running 
coupling in this theory have indicated the presence of a 
strongly-coupled infrared fixed point lfT6l [TTl with mass 
anomalous dimension 7"^ ~ 1 IfTSl . Our simulation results 
are found to be consistent with infrared conformality, up to 
large finite -volume corrections which appear at our light- 
est fermion masses Trif < 0.010. A global fit using the 
framework of mass-deformed conformal perturbation the- 
ory ll9l[T9lg2l yields a best-fit value 7* = 1.10(17) (with 
only statistical error shown). 

The mass-deformed conformal fit carried out here is 
rather simple, for example setting g{fj,) ~ g* and thus ig- 
noring weakening of the gauge coupling at energies near 
the UV cutoff. Our numerical estimates for 7* should 
therefore be regarded with caution. More importantly, we 
cannot rule out the possibility that at some energy be- 
low the induced confinement scale M corresponding to 
~ 0.015, intrinsic confinement and spontaneous chi- 
ral symmetry breaking would set in. The mass-dependence 
of all observables would then be expected to switch con- 
tinuously from power-law scaling to the forms expected in 
chiral perturbation theory. 

In future work, simulations at additional fermion masses 
and on additional volumes will be crucial. If the Nf = 
10 theory is indeed conformal or near-conformal in the 
infrared, all bound-state masses are expected to fol- 
low universal scaling functions in terms of the variable 
,^1/(1+7*)^ lHOlEIllia. Observation of "curve collapse" 
in measurements on several volumes can therefore provide 
a clearer signal that the underlying theory is described by 
an infrared fixed point, at least over the range of energy 
scales considered. 

It is instructive to compare the situation at Nf = 10 
to the current state of knowledge for the Nf = 12 the- 
ory. Although the spectrum-fitting methods used here can- 
not distinguish a truly conformal theory from one with a 
small dynamical breaking scale /x ^ 1 /L, they do lead to 
7* > 0.8, matching expectations for a theory near the edge 
of the conformal window. For Nf = 12 the same statement 
applies, except that mass-deformed fits indicate 
7* ~ 0.4 <^ 1. Assuming that 7* increases monoton- 
ically as Nf approaches Nf and that the conformal win- 
dow closes near 7* = 1 iflSVETl . a value of 7* ~ 0.4 
indicates that the theory is inside and not particularly close 
to the edge of the window. Furthermore, several groups 
have carried out direct studies of the renormalization-group 
flow, through running coupling studies or similar tech- 
niques H |5j |8j [TTHT31 . The observation of backwards 
renormalization-group flow by some of these groups gives 
additional evidence that the Nf = 12 theory is in the con- 
formal window. 
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Appendix: data tables 

To appear in supplementary journal material. 



0.010 
0.015 
0.020 
0.025 
0.030 



Mp 
0.1954(44) 
0.2023(77) 
0.2252(23) 
0.2482(19) 
0.2762(16) 



My 

0.2442(30) 
0.2443(40) 
0.2660(40) 
0.2855(19) 
0.3185(24) 



Ma 
0.2620(46) 
0.2836(42) 
0.3257(23) 
0.3487(44) 
0.4064(37) 



Mn 
0.3804(59) 
0.3652(79) 
0.3933(114) 
0.4187(55) 
0.4703(57) 



Mn* 
0.4277(69) 
0.4319(18) 
0.4865(80) 
0.5075(107) 
0.6043(115) 



nii_ 



TfXres 



Fa 



(V'VXxiO" 



0.010 
0.015 
0.020 
0.025 
0.030 



0.001737(3) 
0.001765(4) 

0.001812(2) 
0.001855(4) 
0.001883(4) 



0.01720(30) 
0.02289(41) 

0.02852(63) 
0.03204(40) 
0.03772(56) 



0.03611(71) 
0.03551(73) 

0.03728(88) 
0.03811(48) 
0.04193(54) 



0.03479(49) 
0.03111(80) 

0.03188(112) 
0.02838(54) 
0.03137(69) 



1.37790(44) 
1.96148(54) 

2.54567(85) 
3.12672(47) 
3.70510(73) 



TABLE II. Ordered- start value for all Nf = 10 observables. 



rrif 


Mp 


Mv 


Ma 


Mn 


Mn* 


0.010 


0.1829(30) 


0.2418(51) 


0.2686(60) 


0.3637(58) 


0.4107(88) 


0.015 


0.1889(40) 


0.2376(59) 


0.3129(34) 


0.3630(90) 


0.4508(166) 


0.020 


0.2464(28) 


0.2902(24) 


0.3659(24) 


0.4297(48) 


0.6118(181) 


0.025 


0.2673(12) 


0.3209(20) 


0.4190(50) 


0.4774(47) 


0.8392(706) 


0.030 


0.2844(12) 


0.3379(21) 


0.4372(64) 


0.5053(49) 


0.6710(174) 


mf 


TYlres 


Fp 


Fv 


Fa 


(^,A)(xio-2) 


0.010 


0.001740(4) 0.01981(123) 0.03755(83) 0.03608(139) 


1.37796(42) 


0.015 


0.001759(5) 


0.02704(45) 


0.03611(91) 


0.03193(71) 


1.96729(61) 


0.020 


0.001790(4) 


0.03497(98) 


0.04141(59) 


0.03341(29) 


2.53830(53) 


0.025 


0.001820(3) 


0.03794(71) 


0.04634(91) 0.03628(122) 


3.11588(104) 


0.030 


0.001869(2) 


0.03878(71) 


0.04502(68) 0.03558(142) 


3.69653(82) 



TABLE III. Disordered-start value for all Nf = 10 observables. 



ruf 


Mp 


Mv 


Ma 


Mn 


Mn* 


0.010 


0.1891(62) 


0.2430(92) 


0.2653(150) 


0.372(15) 


0.419(60) 


0.015 


0.1960(80) 


0.2410(95) 


0.2982(162) 


0.364(16) 


0.441(65) 


0.020 


0.2359(66) 


0.2781(99) 


0.3458(184) 


0.411(18) 


0.549(79) 


0.025 


0.2578(69) 


0.3032(104) 0.3838(208) 


0.448(18) 


0.673(108) 


0.030 


0.2803(74) 


0.3282(113) 0.4218(229) 


0.488(19) 


0.638(92) 


mf 




Fp 


Fv 


Fa 


(^V^)(xl0-2) 


0.010 


0.001738(10) 


0.0185(14) 


0.0368(22) 


0.0354(27) 


1.3779(16) 


0.015 


0.001762(11) 


0.0250(15) 


0.0358(22) 


0.0315(23) 


1.9644(21) 


0.020 


0.001801(10) 


0.0317(20) 


0.0393(23) 


0.0326(24) 


2.5420(27) 


0.025 


0.001837(11) 


0.0350(21) 


0.0422(25) 


0.0323(25) 


3.1213(32) 


0.030 


0.001876(11) 


0.0383(23) 


0.0435(25) 


0.0335(26) 


3.7001(38) 



TABLE IV. Combined values for all Nf = 10 observables, using 
the method described in the text. 



